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ABSTRACT 

We present a scenario for the production of dynamically cold rings around the Milky Way via a high 
eccentricity, flyby encounter. These initial conditions are more cosmologically motivated than those 
considered in previous works. We find that the encounters we examine generically produce a series 
of nearly dynamically cold ring-like features on low-eccentricity orbits that persist over timescales 
of ~ 2 — 4 Gyr via the tidal response of the primary galaxy to the close passage of the satellite. 
Moreover, they are both qualitatively and quantitatively similar to the distribution, kinematics, and 
stellar population of the Monoceros ring. Therefore, we find that a high eccentricity flyby by a 
satellite galaxy represents a cosmologically appealing scenario for forming kinematically distinct ring- 
like features around the Milky Way. 

Subject headings: galaxy: kinematics and dynamics, galaxy: structure, galaxies: interactions, galaxies: 
structure, methods: n-body simulations 



1. INTRODUCTION 

We live in a hierarchical unive rse, in which merg- 
ers are a frequent occurrence (lLacev fe Cold li~99l 
iSomerville fc Kolattl H99l iSomerville et alJ 1200(1 ). In 
this scenario, galaxies like the Milky Way (MW) build up 
much of their mass by accreting smaller satellite galax- 
ies. These minor mergers are connected to a variety of 
observable stellar structural and kinematic phenomena in 
the MW and exter nal galaxies, including: the build-up 
of ste llar halos ("e.g. lBullock fc Johnstonll2005HBell et al l 
20071). so-called "anti t runcated" (lErwin et al.l 120051 : 
Poh len fc Trujilld 120061: lYounger et all I2007D and 



ex- 



tended" outer disks (llbata et al.ll2005t IPenarrubia et al.l 
120061 : llbata et al.l 120071) . and the dynamical heating of 
the st ellar disk (e.g iToth fc Ostrikerlll992HQuinn et all 



fT99l IWalker "et alJfl996l IKazantzidis et al.ll2007l) . 

In our own galaxv JNewberg et al.1 ( 20021 ) recently iden- 
tified a coherent ring-like structure in Monoceros (MRi), 
using data from the Sloan Digital Sky Survey (SDSS: 
I York et al.ll2000D . Since this initial detection, the MRi 
has been identified in the infrar ed from the Two Mi- 
cro n All-Sky Survey (2MASS: ISkrutskie et~aI1 [2006h 
bv iRocha-Pinto et alJ (|2003l) and IMartin et all (2 004). 
reanalyzed using upd ated SDSS data (|Grillmairi 120061 : 
iBelokurov et al.ll2006f ). and has been followed-up photo- 
metrically by seve r al op tical surveys (jlbata et al.l 120031 : 
iConn et alJl2005bL 12007ft . The MRi subtends £ 100° 
in galactic longitude and lies at a galactocentric radial 
distance of between 15 and 20 kpc. Spectroscopic stud- 
ies have shown that it is kinematically distinct from the 
disk - dynamically cold with a low eccentricity orbit - 
and is composed primarily of low-metallicity stars with 
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1.6 <; [Fe/H] <> - 0.4 (iCrane et alJ 



120031 : IConn et al.ll2005al IMartin et al 



2003- lYanny et all 
20051 12006ft ~ 



Numerical modeling has suggested that similar struc- 
tures could be formed via tidal disrup tion of a low-mass 
satellite on a nearly coplanar orbit (Helmi et al. [20031: 



Meza, et all 120051: IMartin et all 120051 : IPenarrubia et al I 
20051 120071) . The most successful of these models ar- 



gue for a very low mass companion with M s /Mmw ~ 
1/1000 on a nea rly ci r cular (e = 0.1) progr ade or- 
bit ([Martin et al.l 120051 : IPenarrubia et al.l 1 120051. 120071) . 
However, recent cosmolog ical simulations (|Bensonl 120051 : 
iKhochfar fc Burkertll2006ft have found that the orbits of 
satellite galaxies are likely be highly eccentric. The dy- 
namical friction timescale for such interactions is very 
long and as a result it is unlikely that such a satel- 
lite could have circularized in time to form the MRi 
(|Besla et al.l l2007t iBovlan-Kolchin et al.l l2008f ) . There- 
fore, while they are successful at reproducing many of 
the characteristics of the MRi, these models may have 
cosmologically unappealing initial conditions. 

With this in mind, we propose an alternative mech- 
anism for forming dynamically cold ring-like structures 
around the MW: a flyby encounter with a small satel- 
lite on a hig h eccentricity orbit . This scenario was first 
suggested bv IKazantzidis et al.l (|2007l ), but here we in- 
vestigate the triggering event in detail, and discuss its 
relevance to the MRi. 

2. SIMULATIONS 

The simulations pre sented in this study were per- 
formed with GADGET2 (|Springell2005ft . anN-Body/SPH 
(Smooth Particle Hydrodynamics ) code using the en- 
tropy conserving formalism of Sprin ge! fc Hernquistl 
(2002). We include radiative cooling and star fo rmation, 
tuned to fit the observ ed local Schmidt Law (| Schmidt! 
119591 : iKennicutH [1998). We also incorporate a sub- 
resoluti on multi-phase model of th e interstellar medium 
(ISM) (|Springel fc Hernquistl l2003ft - softened {q Eo s = 
0.25) such that the mass-weighted ISM temperature is 
~ 10 4 ' 5 K - and sink particles representing supermassive 
black holes that accrete gas and release is otropic thermal 
feedback to self-regulate their growth (Springel et al. 
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Fig. 1. — The evolution of the projected stellar mass density, colored according to a logarithmic scale. The panels are 70 kpc on a side, 
and the simulation time is printed in the upper left hand corner in units of Gyr. The color-bar indicates the projected K-band surface 
brightness in units of mag arcsec -2 , assuming a constant M/L ratio of 2 in solar units. 

Star formation induced by the merger is co ncentrated 
in the nucleus of the primary galaxy (see e.g.JHernciuist 
[i989|lMihos fc HernauistJll994[l996l:lHernauist & Mihos 
1 1995! ) . making these rings primarily disk stars. They are 
furthermore dynamically cold and kinematically distinct 
from typical disk stars (see Figure[3|), and slowly disperse 
on a timescale of ~ 2 — 4 Gyr (or several rotation periods) 
owin g to phase mixing of the collisionlcss stellar parti- 
cles (jBinnev fc Tremaindll987l ). Because these rings are 
generated via gravitational interactions, they are largely 
insensitive to the gas content or structural parameters of 
the satellite galaxy. And because the flyby is effectively 
an impulse interaction, they will still be formed even if 
the satellite is disrupted during the encounter. 



I2005blh 

The progenitor galaxy mod els were constructed fol- 
lowing ISpringel et alj (|2005af ). to which we refer the 
reader for details. The primary galaxy (total mass of 
M200 = 10 12 /i _1 M Q ) is analogous to the MW with a ba- 
ronic mass fractio n of mi, = 0.05. I t was realized with 10 6 
halo particles in a lHernquisll (|1990f ) profile with a concen- 
tration of c = 9 as m otivated by cosmological simulations 
([Bullock et al.ll200l . and 4 x 10 5 stellar disk (80% of 
the baryonic mass) and 2 x 10 5 gas particles (20% of the 
baryonic mass). The satellite galaxy had a total mass of 
M 2 oo = 5 x lO 10 /^ 1 Mq (Mpg/Msg ~ 20) and an iden- 
tical baryonic mas s fraction. It wa s realized with 10 5 
halo particles in a iHernquistl (|1990) profile with a con- 
centration of c = 18, again motivated by cosmological 
simulations, and 2 x 10 4 stellar disk (50% of the bary- 
onic mass) and 4 x 10 4 gas particles (50% of the baryonic 
mass). They were placed on a parabolic encounter with 
Bp = 5/i _1 kpc perigalactic radius (~ 1 scale length), 
consistent with the results of cosmologi cal simulations 
(|Bensonll200"l iKhochfar fc BurkertJl2006f) . 

3. DISCUSSION 

The interaction is summarized in Figures [T] and 
[21 After the initial close passage, resonances be- 
tween the orbital frequency of the satellite and pri- 
mary galaxy disk particles ( "sta rs" ) excite coplanar tidal 
arms ( Toomre fc Toomrel [1972). These features then 
wrap around the primary galaxy as it continues to re- 
volve, forming a set of concentric ring-like features. 



4. RELEVANCE TO THE MONOCEROS RING 

The ring features that are produced by the flyby inter- 
action ar e similar in several wa ys to the MRi first iden- 
tified by iNewberg et al.l (|2002t ). They provide a good 
match to the kinematics and location of the ring, and 
are roughly consistent with metallically measurements of 
MRi stars. And, cosmological simulations indicate that 
such interactions are common for MW-sized halos over 
the t imescale for ring form ation (i.e., within the past ~ 4 
Gvr: iStewart et aT1l2007l ). Therefore, while we do not 
claim to have modeled the MRi in detail, based on these 
similarities we propose a flyby interaction as a possible 
formation scenario. 

The MRi forms a coherent structure over ~ 
100° in galactic longitude from approximately 15-20 
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Fig. 2. — Three different projections (X— Y, Z-Y, and X-Z) of the K-band surface brightness (left; again assuming a constant M/L ratio 
of 2 in solar units) and gas surface mass density (right) at t fa 3 Gyr. Panels are 50 kpc on a side in the X and Y directions, and 20 kpc 
in the Z direction. The solid line corresponds to one of the ring structures as roughly the same galactocentric location (R ga c = 19 kpc) as 
the MRi. 
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kpc from the galactic center dNewberg et al.l l2002t 
Ibata et alj|2003l; |Rocha-Pinto et all I2003I: iMartin et all 



20041 IConn et al.ll2005bll2007f ). Kinematicallv. it is dis- 



tinguished from MW disk stars by its low radial velocity 
dispersion (~ 25 km s" 1 ); the MRi is dynamically cold 
(ICrane et al.||2003| : lYannv et al.ll2003t IConn et al.ll2005al : 
IMartin et al.l l2005. 2006:). Proper motion measurements 8 



Fig. 3. — Phase space diagram for the stellar particles at t f» 3 Gyr, showing the orbital frequency (O) as a fraction of that expected from 
a circular orbit at that radius (O fe ) as a function of galactocentric radius (i? gac ). The left panel shows the density of points in gray-scale, 
with linear contours overlaid to guide the eye and a dashed line at the location indicated in Figure f3] The right hand panel is color coded 
by the radial velocity dispersion (<t(vji)) in each cell, with the same contours overlaid. 

- with a similarly low radial velocity dispersion (see Fig- 
ure [3]) ■ This ring also extends —4 <^ z <J 4 kpc above 
and below the galactic pl ane, which i s cons istent with 
detections such as those of IConn et al.l ()2007l ). Although 
it is unclear what fraction of the ring would be identified 
in observations, it contains <J 1% ( < 5 x 10 8 M Q ) of the 
total stellar mass of the disk, which is consistent with 
the mass estimates of lYannv et al.l ()2003l ). 

A disk origin for MRi stars is also broadly consistent 
with their observed stellar population . Several spectro- 
scopic studies (e.g.. lYanny et al.ll2003l ) have found mean 
metallicities of [Fe/H] = — 1.6 ± 0.3, while ICrane et al.l 
(2003) find a higher mean metallicity of [Fe/H] « -0.4 
using a different tracer population. This suggests a pri- 
marily metal-poor stellar population with a spread in 
metallicities, and possibly multiple epochs of star forma- 
tion. Our modeling indicates that the MRi may have 
formed from outer disk stars moved outwards by tidal 
interactions with the satellite galaxy, whic h is consis- 
tent with this metal-poor stellar population (|Luck et al.l 



also sugg est that the ring stars are in low-eccen tricity or- 
bits (e.g.. ICrane et al.ll2003HYannv et al.ll2003fh We find 
that our simulations produce a ring in the right galac- 
tocentric radial range that is approximately circularly 
supported - 0.9 ^ il/ilk ^ 1-1, where f2 is the orbital 
frequency and Qk is that expected from circular motion 



8 Crane et al. (2003) measure an overall circular velocity of 
~ 100 km s — 1 . However, we note (as in [Pcnarrubia et afl 120051) 
that the accuracy of such measurements in physical units is lim- 
ited to Av perp ^ 4.74RqA/x where Rq as 8 kpc is the galacto- 
centric radius of the sun and A/t ~ 3 — 4 mas yr~ 1 is the typical 
proper motion measurement error. We have therefore chosen to 
interpret these measurements as being generally consistent with 
low— eccentricity orbits. 
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Fig. 4.— Three projections (X-Y, Z-Y, and X-Z) of the K-band 
surface brightness (again assuming a constant M/L ratio of 2 in 
solar units) at the same simulation time as Figures [2] and [3] for a 
Msg/Mpq ~ 1/100 interaction with identical orbital parameters. 



120061 lYong et all [2006). Furthermore, the tidal inter- 
action moves a significant supply of cold gas into the 
same ring structures (see Figure 0) , which provides the 
raw material for subsequent epochs of star formation. 
While the density-dependent prescription in our simu- 
lations does not accurately capture some modes of star 
formation owing to the limitations of the SPH method, 
prescription s that include an approximate treatment of 
shocks (e.g.. lBarnesll2004[ ) may show these multiple star 
formation episodes. 

While ring features in our simulations capture many 
of the properties of the MRi, there are some observa- 
tions which potentially conflict with our proposal. In 
particular, the satellite in our simulations is roughly the 
same mass as the LMC, and by ~ 2 Gyr after the in- 
teraction would be at a galactocentric distance of ~ 250 
kpc. Such a massive object is well wit hin the detec- 
tion t hreshold of current surveys (e.g., IWillman et al.l 
l2002t IKoposov et al.1 [20071 ); Leo I - a far lower mass 
MW companion - was detected at compar able distance 
(ICaputo et al.lll999t iBellazzini et al.l 120041 : IMateo et all 
|2007|) . However, dynamically cold ring structures are 
produced generically in fly by encounters, and similar 
features are present in both a lower-inclination inter- 
action (i — 10°) and for a lower mass satellite galaxy 
(Ms G /Mp G ~ 1/100; Figure [4]). It is thus possible that 
the satellite galaxy that produced the MRi is either (1) 
hidden from view, at low galactic latitude and/or on the 
other side of the galaxy, or (2) it is lower mass with 



1/100 <> Msg/Mpg ^ 1/20. Therefore, it is not unrea- 
sonable to speculate that Leo I, which had at least one 
assage through the MW disk within the past 2-4 Gyr 
Sohn et aTll2007f) . could have excited the MRi we see to- 
day. Moreover, it is likely that this interaction stripped 
a significant fraction of its mass 9 , making current dy- 
nami cal estimates (Ml bo i/Mmw ~ 1/1000: IMateo et all 
|2007f ) a lower limit on its mass at the time of the in- 
teraction. Furthermore, simulations indicate that there 
are likely > several sub-haloes in the MW halo with 
1/100 <; Msh/ Mpq <s 1/20 on orbits that take then 
through the disk riMoore et al.lll999t iGiocoli et"alll2007t 
iKazantzidis et ai1 l2007). Since gala xy formation may be 



inefficient in low mass haloes (e.g.. iHaiman etall[l997t 
iBarkana fc Loeblll999l iKravtsov et al.ll2004h . thev may 
be significantly underluminous and therefore difficult to 
identify. 

Finally, there ha v e been suggestions in t he literature 
([Crane et al.l 120031 : iFrinchabov et alj I2004D that some 
globular clusters (GCs) may be associated with the MRi. 
This is also potentially in conflict with our modeling. 
However, we note that their physical association with 
the MRi is somewhat speculative. Alternatively, there 
are ^s everal nearly coplanar GCs at roughly the same 
radius (|Harrislll996D as progenitor MRi stars in our sim- 
ulations that they could have also been moved into the 
rings by the interaction. 

5. CONCLUSION 

We present a scenario for the production of dynami- 
cally cold rings around a disk galaxy such as the MW 
via a prograde flyby encounter with a satellite galaxy. 
Tidal arms excited during close passage coalesce and 
wrap around the disk of the primary galaxy. These kinds 
of interactions are more cosmologi cally likely than the 
nea rly circular orbits presented b y iMartin et al.l (2005) 
and lPenarrubia et al.1 (|2005l [20071 ). while dynamical fric- 
tion is insufficient to circularize the orbit of such a low 
mass companion. Our modeling shows rings with simi- 
lar spatial distribution and kinematics to the MRi. The 
disk origin of MRi stars is furthermore broadly consis- 
tent with observed stellar populations. Therefore, we 
find that a flyby encounter represents a more cosmolog- 
ically appealing scenario for the production of the MRi 
and other dynamically cold rings around the MW. 
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9 In our simulations, both the 1/20 and 1/100 encounters 
stripped ~ 50% of the satellite mass. While the efficiency of this 
process is sensitive to both the initial structure of the satellite and 
the orbital parameters of the interaction, this suggests that a a 



passage through the disk is likely to remove a significant fraction 
of the satellite's total mass. 
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